The dependence of topological Anderson insulator on the type of disorder 
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This paper details the investigation of the influence of different disorders in two-dimensional 
topological insulator systems. Unlike the phase transitions to topological Anderson insulator induced 
by normal Anderson disorder, a different physical picture arises when bond disorder is considered. 
Using Born approximation theory, an explanation is given as to why bond disorder plays a different 
role in phase transition than does Anderson disorder. By comparing phase diagrams, conductance, 
conductance fluctuations, and the localization length for systems with different types of disorder, 
a consistent conclusion is obtained. The results indicate that a topological Anderson insulator is 
dependent on the type of disorder. These results are important for the doping processes used in 
preparation of topological insulators. 
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I. INTRODUCTION 

With research breakthroughs in HgTe/CdTe quantum 
wells and 3D topological materials^— topological insula- 
tors have attracted much attention in recent years. These 
unique topological properties are responsible for some 
interesting and surprising phenomena. For example, in 
2009, Li, Chu, Jain and ShenS discovered that Anderson 
disorder can lead to a topological phase transitions with 
quantized conductance and named this phase topological 
Anderson insulator (TAI). The phenomenon of TAI was 
also reported by Jiang et al£ in their work. Subsequently, 
the origins of TAI, as well as TAI in other systems, have 
been studied by many groups using a variety of meth- 
ods^-i£ 

Guo et aJjiS found that the TAI phenomenon also ex- 
ists in disordered 3D topological insulator. Recently, 
Xing et aJfii compared the disorder effects in three dif- 
ferent systems where the quantum anomalous Hall effect 
exists. They observed TAI in these three different sys- 
tems and demonstrated that increasing disorder strength 
produces the TAI phenomenon. Furthermore, Yamak- 
age et al~2. reported similar phase transitions in disor- 
dered Zi topological insulators when considering s z non- 
conserving spin-orbit coupling. 

The origins of TAI have been studied by many groups 
as well. After the initial reports of TAI^ Jiang et al£ 
calculated the distribution of local currents in real space 
at various strengthes of Anderson disorder and provided 
an explanation on TAI phase by studying the helical edge 
states in HgTe/CdTe quantum wells. In addition, it has 
been shown by Groth et alM- that the phase transition to 
TAI should ascribe to a negative correction to the topo- 
logical mass because of Anderson disorder. Using the 
effective medium theory and through numerical calcula- 
tions, Groth et al— discussed in detail how Anderson 
disorder renormalizes the topological mass, chemical po- 
tential and finally induces a phase transition. This inter- 



pretation provides a clear physical image of TAI and has 
been accepted generally by most physicists in this field. 

The TAI effect has also been identified in 3D topo- 
logical insulators and the effective medium theory noted 
above is found to be useful when describing the 3D case as 
reported in the paper by Guo et al^ However, by looking 
at the phase diagrams, Prodani^ takes the position in his 
paper that TAI should not be considered a distinct phase 
but should be described as part of the quantum spin-Hall 
phase. In addition, Chen et alM- discussed the TAI phe- 
nomenon from the point of band structures. Finally, by 
calculating the Z2 topological number in systems with 
periodic disordered supercell regimes, Zhang et al*£ ver- 
ified that TAI corresponds to a topologically non-trivial 
phase. 

Thus, it has been well documented that Anderson dis- 
order may induce TAI. However, it is not clear whether 
TAI is certain to be observed experimentally and whether 
all types of disorder can definitely push a phase transi- 
tion to TAI for an anomalous quantum Hall system. It 
is the purpose of this paper to address these questions. 

Besides the on-site Anderson disorder, there exists an- 
other type of disorder, bond disorder, which originates 
from the deformation of the lattice or from some other 
interactions that induce a random hopping term. This 
type of disorder exists widely and cannot be ignored 
when describing a real system. This is especially the 
case when considering a 2D system such as graphene. 
We note that such bond disorders have been studied ex- 
tensively in various systems*^— Following the effective 
medium theory^ it is shown that, unlike Anderson disor- 
der which appears through gq or o z term in Hamiltonian, 
bond disorder, which appears through g x or a y term in 
Hamiltonian, would renormalize the topological mass by 
adding a positive mass correction. Therefore, according 
to the effective medium theory, TAI phenomenon cannot 
arise. To further investigate this issue, we will describe 
two different models. One is the HgTe/CdTe quantum 
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wells and the other is the Haldane model^ These mod- 
els have been chosen because the first model is always 
used to study the TAI phenomena, and in the second 
model the a x (a y ) disorder can be easily introduced by 
the deformation of a honeycomb lattice and the cr x (a y ) 
disorder comes very naturally if there is a random hop- 
ping correction to the nearest hopping term. We have in- 
vestigated the phase diagram, conductance, conductance 
fluctuations and the localization length for these two con- 
crete models. It will be shown that all of the results are 
consistent and that a x (a y ) disorder cannot lead to TAI 
and this type of disorder prohibits the TAI phenomenon 
in some sense. Therefore, in a real system, the presence 
of TAI phenomena may be determined by which type of 
the o x ((jy) and <7o(<7 2 ) disorder is stronger. 

The rest of this paper is organized as follows: In Sec. ITII 
we introduce the two models in the tight-binding repre- 
sentation and derive the formulas of the conductance, the 
renormalized topological mass M and the renormalized 
chemical potential JI. The numerical results are discussed 
in Sec. IIII1 Finally, a brief summary is given in Sec. IIVI 



II. THEORETICAL MODELS 



the lattice Hamiltonian [Eq. (1)] in k-representation 
can be reduced to the continuous Hamiltonian in Ref.6 
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Figure 1: (Color online) A schematic diagram of a infinitely 
long ribbon on the square lattice for the HgTe/CdTe model 
(a) and a infinitely long ribbon in the honeycomb lattice along 
the zigzag direction for the Haldane model (b). 



The first model, which has been studied extensively, 
is the standard HgTe/CdTe quantum wells. The tight- 
binding Hamiltonian for the square lattice sketched in 
Fig. 1(a) has the form^ 



Hi = 



Ed 



Ed 



E s Ui 
Ui Ep 
E„ Ui 







Ui E 



p 



V* 






Vss 

-V. 



Sp 



ri 



\ 



Vsp 

Vpp 



( Vss iV sp 

iv; p Vpp 



\ 






v.. 



(fii+Sx + h.c. 






-iV 



iVcr) Vry, 



Sp 



fi+Sy + h.C. 



(1) 



Here a is the lattice constant and all of the parameters 
A, B, C, D, and M can be controlled experimentally^. 
The topological mass M can be tuned continuously by 
changing the thickness of the HgTe and subsequently 
switches the HgTe/CdTe wells between a topologically 
nontrivial phase and a topologically trivial phase. In 
this model, the individual spin-up Hamiltonian and spin- 
down Hamiltonian in Eq. (1) are time- reversal symmetric 
to each other. Since they are decoupled, we can deal with 
them individually. For simplicity, we shall focus only on 
the spin-up Hamiltonian in the following calculations. 

The Haldane model proposed by Haldane in 1988^ 
considers honeycomb lattice with next-nearest-neighbor 
coupling and a staggered sublattice potential. The 
Hamiltonian can be expressed as^ 
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Here i = (i^, i y ) is the site index, and Sx and Sy are unit 
vectors along the x and y directions, (pi = (a[,Ci,bi,di) T 
represents the four annihilation operators of the electron 
on the site i with the state indices |s,t))|Pa: + ipy,t) > 
\s, |), | - {p x - ip y ), i> respectively. E Sl E p , V ss , 



Vpp, and 



V sp are the five independent parameters that character- 
ize the clean HgTe/CdTe samples. Ui represents ran- 
dom bond disorder, which is uniformly distributed in the 
range [— tt] with the disorder strength U ^2r— Note 
that, in real materials, the disorder strength in same 
sites should be much stronger than that between neigh- 
bor sites and therefore only random bond disorder for the 
same cell is included^ It is clear that near the T point, 



The first term is the onsite energy and e m = ±A for the 
A site (B site), which is shown by a dot (circle) in Fig. 
1 (b) . The second term is the usual nearest neighbor hop- 
ping term. Here, random bond disorder is introduced by 
tmn = tnm = t\ + St r , where St r is uniformly distributed 
in the range [—■§■, ^] with the disorder strength U. Note 
that the random hopping term gives the bond disorder 
a more concrete physical image of the underlying nature 
of the term. The third term is the second neighbor hop- 
ping term with bond dependent phase. Note that v mn 
is different for different hopping directions ^ 11 ' 22 For ex- 
ample, if an electron in A or B site makes a left (right) 
turn to get to the second site, v mn = !(—!)■ 
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By performing Fourier transformations, we can easily 
obtain the low energy effective Hamiltonian for the two 
models given above. Note that the low energy effective 
Hamiltonian for the two models have the same represen- 
tation and can be written as follows^ 

H(k) = H (k) + V = a(k x cr x — ky(Ty) 

+ (M + f3k 2 )a z + (v + jk 2 )a + U(r)a d , (3) 



where the operator k x / y can be represented with the 
momentum operator: k x / y = —ihd x / y . This Hamilto- 
nian is a two-dimensional Dirac Hamiltonian where the 
three Pauli matrices a x , a y , a z and the unit matrix ao 
represent the pseudospins. Note that the pseudospin is 
formed by the s or p orbital (corresponding to the A or 
B sublattice) for the HgTe/CdTe (Haldane) model. The 
scalar potential U(r) accounts for the disorder amplitude 
and ad = a / x denotes Anderson disorder or bond disor- 
der respectively. The parameters a, fx, 7 and /3 for the 
HgTe/CdTe model have the simple form of: 



-A; fi = C; 7 = -D; [3 = -B. 



(4) 



Meanwhile, the low energy effective Hamiltonian for 
the Haldane model can be expanded at the two inequiv- 
alent Dirac points K and K* . Therefore, the parameters 
a, ji, M, 7 and j3 for the Haldane model can be repre- 
sented as: 



SoCo+E^a^+Ei/Cry-l-Eza'z. Obviously, the role of T>(E, k) 
induced by ao or &x disorder must be to make a correc- 
tion to the <jj terms in H(k). The topological mass and 
the chemical potential are then renormalized and take 
the form: 



M = M 



lim i?eS 2 

fe^O 



Ji = E 



lim ReZ - (7) 



To acquire self-energy in numerical calculations, 

the self-consistent Born approximation is applicable and 
the integral equation for T,(E) and can be written as: 

= ^-(a/2Tr) 2 [ d 2 k[a d (E+I - H „ - £)-V d ],(8) 
Jbz 

where ad = a / x denotes the type of disorder. 

Up this point, it becomes apparent that the different 
types of disorder may lead to a different correction to the 
topological mass, the sign of which is critical for classi- 
fying the topological phase of the system. To see this 
clearly, we neglect on the right side of Eq. (8). An 

approximate solution for E(-E) with a closed form is then 
given by: 
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(9) 
(10) 



a = —t\a\ fi = —\/3t2Cos(j>; M = A ± 3V3t2sin<f>, 



9 9\/3 
7 = — t2d 2 cos<f>; P — =F t2a 2 sin(f>, (5) 



where ± or =p corresponds to K or K* respectively. 

It has been explained^ that elastic scattering by onsite 
Anderson disorder (called ao disorder to distinguish it 
from bond disorder being called a x disorder correspond- 
ingly in this paper) causes a state with a definite momen- 
tum to decay exponentially as a function of space and 
time. Therefore, a negative correction to the topologi- 
cal mass M is induced due to the quadratic term in the 
Hamiltonian. The renormalized mass M = M + 6M may, 
in fact, have the opposite sign to the bare mass M. This 
implies that onsite Anderson disorder, erg, may induce 
a phase transition from a normal insulator to the topo- 
logical phase. When this happens, TAI phenomena can 
be observed. However, the question arises as to whether 
this applies to the bond disorder, namely a x disorder. 
We will now qualitatively show this mechanism through 
the following derivation. 

First, we begin with the Hamiltonian in Eq. (3). Fol- 
lowing the same derivation given in the paper by Groth et 
al.r^ the self-energy can be obtained from the equation: 



\E-H a {k)-^{E,k)]- 1 = ([E - H(k)]' 



(6) 



where (. . .) denotes the disorder average. Then, the 
self-energy can be represented in the form: £(£7, k) = 



where =p corresponds to <7q and o x disorder respectively. 
For random bond disorder, the cto term in the self- 
energy T,(E,k) above is obtained by o x o$o x = Co; f° r 
normal Anderson disorder, it may be calculated using 
"WoCo = °o- Thus, the a term does not produce any 
difference for the two types of disorder, and the Fermi 
energy, which corresponds to the ao term, has the same 
renormalization for normal Anderson disorder, ctq and 
bond disorder, a x . However, the a z term in the self- 
energy Ti{E,k) is changed differently by Anderson dis- 
order and bond disorder. Because aoa z ao = &z and 
a x a z a x = — er z , the topological mass, which corresponds 
to the a z term, is renormalized along opposite directions 
by Anderson disorder ao and bond disorder a x . This can 
be seen clearly in Eq. (9). Hence, TAI phenomenon in 
a a x disordered system may manifest different features 
from that in a ao disordered system. 

Although the above formula are derived in the 
HgTe/CdTe quantum wells, they should be also appli- 
cable to the Haldane model since the low energy effective 
Hamiltonian for the two models have the same represen- 
tation as BHZ models In addition, it should point out 
that the disorder-induced mass inversion always corre- 
sponds to a topological phase transition in HgTe/CdTe 
quantum wells, however it is not true for the Haldane 
model if the topological mass only change its sign at K 
or K* point. That is because the topological property of 
the Haldane model is determined by the relative sign of 
the topologically effective mass at the two Dirac points 
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K and K*. Namely, if the topological masses have oppo- 
site sign at K and K* points, the system is topologically 
non-trivial. Otherwise, the system is topologically triv- 
ial. In the following, we will see clearly that the disorder, 
regardless of Anderson disorder or bond disorder, has the 
same effects on the two models. 

From the models in Eqs. (1) and (2), it is easy to de- 
scribe a nanoribbon geometry as shown in Fig. 1. Here, 
only the nanoribbon with a zigzag edge is studied for the 
Haldane model. Using the Landauer-Biittiker formula, 
the linear conductance at zero temperature and low bias 
voltage can be represented as> & 23 ' 24 

G LR = ^-T=^-Tr[T L G r T R G% (11) 
h h 

where T = Tr[T iGT R G a ] is the transmission coef- 
ficient from the left lead (source) to the right lead 
(drain), T L/R = i(£ r L/R - with being the 

retarded/advanced self energy, respectively. To present 
a more concrete picture, we will provide some numerical 
calculations for the two models and follow with discus- 
sions on the implications in section III. 

III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. Band structures for two models 

In this section, the conductance, conductance fluctu- 
ations and localization length are studied numerically. 
The size of the central region is denoted by integers N 
and W, which represent the length and width respec- 
tively. For example, in the schematic diagram shown in 
Fig. 1(a), the length of the central region (red) is given 
by L x — Nxa with N = 9, and the width by L y = W x a 
with W = 5. In Fig. 1(b), the length of the central re- 
gion (red) is given by L x = N x y/3a with N = 5, and 
the width by L y = W x \fia with W — 3. Here, a is 
the square lattice constant a = 5nm for the HgTe/CdTe 
model and represents the nearest neighboring distance of 
a = 0.142nm for the Haldane model. 

We have plotted the band structures of the two mod- 
els in Fig. 2. Figs. 2(a) and 2(b) are plotted for the 
HgTe/CdTe model with parameters M = lmeV 2(a) 
and M = -lOmeV 2(b); Figs. 2(c) and 2(d) are for 
the Haldane model with A = — 1.6i 2(c), and A = 
— 1.77i 2(d). In the numerical calculations, other sample- 
specific parameters are fixed to be A = 364.5meV-nm, 
B = — 686meV-nm 2 , C = 0, D = — 512meV-nm 2 for the 
HgTe/CdTe model. For the Haldane model, t is set as 
the energy unit (t = 1), and other parameters are set 
with values of t\ = —t, t 2 = — 0M, <j) = — 0.2357T. It 
can be clearly seen that for the topologically nontriv- 
ial phase of two models, gapless states traverse the bulk 
gap in Figs. 2(b) and 2(d). When changing the topo- 
logical mass from the topologically nontrivial phase to 




ka/rc ka/n 



Figure 2: (Color online) The band structures for the 
HgTe/CdTe and Haldane models are shown in (a), (b) and 
(c), (d). (a) and (b) are plotted with different topological 
masses of M = lmeV and M = — lOmeV; (c) and (d) are 
plotted with different staggered terms of A = — 1.77t and 
A = — 1.60i. The red (blue) lines represent the bulk states 
(edge states). 



the topologically trivial phase, gapless states disappear 
from the bulk gap for 2(a) and 2(c). If the stripe ge- 
ometry in Figs. 1(a) and 1(b) has a periodic boundary 
in the y direction, namely transforming the system into 
a cylindrical geometry, the gapless edge states disappear 
in 2(b) and 2(d), but no visible changes can be observed 
in 2(a) and 2(c). For simplicity, the band structures for 
the cylindrical geometry are not shown here. This means 
that the two systems are both 2D topologically trivial 
insulators for the parameters in Figs. 2(a) and 2(c) but 
are topologically nontrivial insulators for the parameters 
in Figs. 2(b) and 2(d). 



B. The conductance and its fluctuations 

In Figs. (3) and (4), we investigate the conductance and 
its fluctuations versus strength of the two types of disor- 
der in the models. In all of these calculations data are 
averaged for 500 random configurations of disorder. Note 
that for Anderson disorder, the results show the same 
physics as that described in previous paper o 8 ' 9 ' 13 for the 
HgTe/CdTe model [see Figs. 3(a) and 3(c)]. The conduc- 
tance and its fluctuations are specifically studied for the 
Haldane model with Anderson disorder [see Figs. 4(a) 
and 4(c)] and the same quantum conductance plateau 
and zero conductance fluctuation can be seen within a 
definite range of Anderson disorder. Thus, for Anderson 
disorder, we find the following properties: First, the con- 
ductance decays and its fluctuations initially increases 
gradually when increasing strength of Anderson disor- 
der. Second, at moderate Anderson disorder strength, 
the conductance stops deceasing and falls to a quantum 
plateau (e 2 //i, when considering only spin up or spin 
down case). In addition, the conductance fluctuation re- 
duces to zero [see Fig. 3(c)] or a very small value [see Fig. 
4(c)] over the corresponding range of Anderson disorder 
strength^ Third, as the strength of Anderson disorder 
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continues to increase, the conductance and its fluctua- 
tions both decrease gradually to zero. The system thus 
finally transforms to an Anderson insulator when the dis- 
order is sufficiently strong. 
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Figure 3: (Color online) The conductance G (a), (b) and 
conductance fluctuations SG (c), (d) vs disorder strength U 
for the HgTe/CdTe model. Different curves correspond differ- 
ent Fermi energy E or different topological masses. The size 
of the central region is set at L x = 200a, L y — 80a. Other 
parameters (A, B, C and D) are given at the beginning of 
Sec. IIIII Here, ctq and a x disorders correspond to Anderson 
disorder and bond disorder respectively. 



4(b) and 4(d). Therefore, it can be concluded that bond 
disorder cannot induce a phase transition to TAI for the 
two models as discussed above in Sec. II. 



S. 0.2- 
O 



1 -.1 
' 1\ 
1 

1 \ 
' \ 

1 \ 

V 

* V 


A=-1.60t, E=1.10t 

A=-1.60t, E=0.75t 

A=-1.60t, E=0.50t 
A=-1.77t, E=0.60t 


t • 
1 J 

'i 
I - 
1 


(b) 

V*. G disorder 

V 


CJ disorder 




1 

I 


— . 




1 \ i / / \j 

i \ 'J >■■ \i 

J , .*■ CS disorder 
\ J ■ 


(d) 

z*- 

i ■* 

i 

| 

* J 


(7^ disorder 

1 i\ - 

! »\ 

i 


1 1 2 3 4 5 6 7 
U(t) 


1 2 3 4 5 6 7 
U(t) 



Figure 4: (Color online) The conductance G (a), (b) and 
conductance fluctuations 8G (c), (d) vs disorder strength U 
for Haldane model. Different curves correspond to different 
Fermi energies E or staggered terms A. The size of the central 
region is set to be L x = 200 x \/3a, L y = 92 x 3a. Other pa- 
rameters (ti, t2, and (f>) are given at the beginning of Sec. IIIII 
Here, ao and a x correspond to Anderson disorder and bond 
disorder respectively. 



In general, the quantum conductance plateau and zero 
conductance fluctuations always imply a new phase or a 
novel phenomenon. From the results about the quantum 
conductance, TAI was initially found three years ago4 As 
mentioned in Introduction, Groth et al— showed how 
Anderson disorder induce a phase transition. Namely, 
Anderson disorder will add a negative correction to the 
topological mass. When the topological mass changes 
its sign at strong Anderson disorder strength, a phase 
transition is triggered from the topologically trivial phase 
to the topologically nontrivial phase. 

However, due to the extensive existence of bond dis- 
order in real materials, it is necessary and important to 
study what happens when considering bond disorder. In 
Figs. 3(b), 3(d), 4(b), and 4(d), the conductance and 
its fluctuations are investigated for the two models with 
bond disorder. Contrary to the case of Anderson disor- 
der, TAI is not observed in the two models when chang- 
ing the strength of bond disorder. For example, for the 
HgTe/CdTe model in Figs. 3(b) and 3(d), a series of 
topological masses M and Fermi energies E are chosen. 
In each case, the conductance gradually falls to zero and 
the quantum plateau shown in Fig. 3(a) disappears com- 
pletely. Moreover, the conductance fluctuations in Fig. 
3(d) only shows a peak and then falls to zero finally be- 
cause of the Anderson localization. This description is 
also qualitatively true for the Haldane model in Figs. 



C. The conductance phase diagrams for two types 
of disorder 

To see this more clearly, the conductance phase di- 
agrams are compared for Anderson disorder and bond 
disorder in the HgTe/CdTe model. The results, in Figs. 
5(a), 5(c), 6(a), and 6(c), are shown for normal Ander- 
son disorder; Figs. 5(b), 5(d), 6(b), and 6(d) correspond 
to bond disorder. It can be seen from Figs. 6(a), and 
6(c) that at moderate Anderson disorder and Fermi en- 
ergy, a clear TAI phase (green region) is present and that 
the TAI phase in Fig. 6(c) must correspond to a nega- 
tive renormalized topological mass M (blue region). The 
Anderson disorder renormalizes the topological mass M 
along the negative direction and Fermi energy along the 
positive direction and therefore induced a phase transi- 
tion to TAI. The results in Figs. 5(a), 5(c), 6(a), and 6(c) 
are similar to those described in previous papersj 8 i 9 i 13 In 
summary, normal Anderson disorder can localize the bulk 
states. For a original system of a topological nontrivial 
phase [Fig. 5(a)], the edge state is more robust against 
Anderson disorder than are bulk states and thus leads 
to a quantum conductance region (green region). For a 
original system with a topologically trivial phase [Fig. 
6(a)], a phase transition from a topologically trivial insu- 
lator to a topological insulator can occur as the strength 
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of Anderson disorder increases, and thus result in a TAI 
quantum conductance region (green region). For Ander- 
son disorder, the topological mass is always renormalized 
by adding a negative correction. A more detailed inter- 
pretation can be found in several previous papersi 8 ' 9 ' 13 
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Figure 5: (Color online) The conductance G (a), (b) and 
topological mass M (c), (d) vs disorder strength U for the 
HgTe/CdTe model with M = -lOmeV. The size of the central 
region is set to be L x — 200a, L y = 80a. Other parameters 
(A, B, C, and D) are given at the beginning of Sec. Mil Here, 
ao and <j x correspond to Anderson disorder and bond disorder 
respectively. 

In the following, we are more concerned about the in- 
fluence of bond disorder on this phenomenon. It can be 
seen from Fig. 5 that at M = — lOmeV, the quantum 
conductance region with a value of e 2 /h (green region) 
above the band gap in Fig. 5(a) disappears completely 
as shown in Fig. 5(b). Namely, the quantum conduc- 
tance region where the value of the conductance is e 2 /h 
(green region) decays to a very small region, and only 
exists within the band gap and at a small strength of 
bond disorder [Fig. 5(b)]. The renormalized mass M, 
which remains negative for the whole phase diagram of 
Fig. 5(c) (blue region), alters its sign between the range 
of 100 meV < U < 150meV in Fig. 5(d). Thus, there are 
very different effects on the mass M for the two types of 
disorder. 

It may be surprising that a larger phase region of bulk 
states (red region) appears in Fig. 5(b) than in Fig. 5(a). 
This implies a greater difficulty in localizing the bulk 
states for the case of bond disorder compared to the case 
of Anderson disorder. This result is understandable if 
the different effects of Anderson disorder and bond dis- 
order are considered. Anderson disorder renormalizes the 
topological mass M along positive direction. That is, the 
effective band gap became larger and thus bulk states are 
shifted upward for the Anderson disorder case. On the 
contrary, the opposite effect appeared for the bond dis- 
order. Because states in the center of the bulk band are 
more difficult to localize than states near the edge of bulk 



band^S a larger phase region of bulk states is observed 
in Fig. 5(b) for bond disorder than that in Fig. 5(a) 
for Anderson disorder. Note that the conclusions above 
are based on the same renormalization of the chemical 
potential ~p for both disorder types as given in Eq. (10). 




U(meV) U(meV) 

Figure 6: (Color online) The conductance G, (a), (b) and 
renormalized topological mass M (c), (d) vs disorder strength 
U and Fermi level E for the HgTe/CdTe model with M = 
lmeV. The size of the central region is set to be L x = 200a, 
L y = 80a. Other parameters (A, B, C, and D) are given at 
the beginning of Sec. Mil Here, <7o and a x correspond to the 
Anderson disorder and bond disorder respectively. 

When the system turns to a topologically trivial phase 
with AI = lmeV as in Fig. 6, there is no indication 
of the TAI phenomenon in the conductance phase dia- 
gram, Fig. 6(b). That is, the conductance region with 
the value of e 2 /h disappears completely for bond disor- 
der. Meanwhile, in the phase diagram with a topological 
mass M = lmeV as shown in Fig. 6(d), the topological 
phase (blue region) gives way to the topologically triv- 
ial phase (red region) and the whole diagram is governed 
by positive renormalized topological mass (red region). 
Therefore, a clear conclusion can be drawn from these 
numerical results that normal Anderson disorder gives 
rise to the TAI phenomenon and bond disorder destroys 
it. In other words, whether the TAI phenomenon can be 
observed depends on a competition between Anderson 
disorder and bond disorder in a system because Ander- 
son disorder and bond disorder renormalize the topolog- 
ical mass along different directions and thus affect TAI 
phenomena differently. 

The above conclusions are also applicable to the Hal- 
dane model. When considering Anderson disorder in Fig. 
7(a) with the topologically nontrivial phase and 7(c) with 
the topologically trivial phase, the conductance phase di- 
agrams show a clear region with quantize conductance Go 
(green region), which is the hallmark of TAI in Fig. 7(c). 
However, for bond disorder, the region of quantized con- 
ductance becomes indistinguishable in Fig. 7(b) for the 
topologically nontrivial phase and in 7(d) for the topo- 
logically trivial phase. This indicates that Anderson dis- 
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Figure 7: (Color online) The conductance G vs. disorder 
strength U and Fermi level E for the Haldane model, (a) and 
(b) are plotted with A = — 1.6t in a topologically nontrivial 
region; (c) and (d) with A = — 1.77f in a topologically trivial 
region. The size of the central region is set to be L x = 200 x 
y/3a, L y = 92 x 3a. Other parameters are chosen with <j) = 
— 0.235-7T, t% = — 0.5i. Here, ao and a x correspond to the 
Anderson disorder and bond disorder respectively. 



order can drive a phase transition to TAI but that bond 
disorder can not. Here, we do not give a topological mass 
diagram similar to Fig. 5(c) or 5(d) for the HgTe/CdTe 
model because the Haldane model is based on the honey- 
comb lattice which means that the topological properties 
of this model are determined by the signs of the effective 
masses at the two Dirac points, as introduced in Sec. [TTJ 
When the signs of the effective masses at two Dirac points 
are opposite, the band at one Dirac point would be in- 
verted from that at another Dirac point, and the system 
is topologically non-trivial. Conversely, when the sign of 
the effective mass at two Dirac points is the same, the 
system corresponds to a topologically trivial phase. In 
addition, the parameter /3 in Eq. (5) has different ex- 
pressions at the two the Dirac points. Consequently, it is 
difficult to present clearly the variation of the topological 
properties as a function of the disorder strength using a 
simple phase diagram of the topological mass. 



D. The localization lengths for two types of 
disorder 

In order to highlight the dependence of TAI on the type 
of disorder, we also plot the localization lengt h 13 ! 26 " — 
X/L y of a 2D ribbon with width L y = W x a for the 
square lattice [Figs. 8(a) and 8(c)] and L y = W x 3a for 
the honeycomb lattice [Figs. 8(b) and 8(d)]. Note that 
the localization length is obtained from the relationship 13 
A = 21imi x ^ 00 L x (ln(G/Go)) _1 by increasing the length 
L x of the system at fixed width L y , namely, L x 3> L y . 
In addition, another calculating method, transfer-matrix 



methodj2S~— is also adopted to guarantee the correctness 
of the localization length. Here we have considered both 
cylindrical geometry and ribbon geometry. 

For the HgTe/CdTe model, the localization length 
X/Ly shown in in Fig. 8(a) first deceases when increasing 
disorder strength for Anderson disorder. With further 
increasing the strength of Anderson disorder, a peak ap- 
pears for the localization length X/L y and then decays 
gradually for both geometries, which can be seen clearly 
in Fig. 8(a). The difference in localization length X/L y 
between the two geometries is that peaks have distinct 
heights and are located at the different strengths of An- 
derson disorder. It can be deduced that because of ex- 
isting edge states, a special metal phase emerges for the 
ribbon geometry near to the position of the peak. How- 
ever, when considering bond disorder, the peak of the 
localization length X/L y vanishes completely as shown 
in Fig. 8(c). This provide a further evidence that bond 
disorder does not induce a phase transition to TAI. As 
is shown in Fig. 8(b) and 8(d) for the Haldane model, 
there is no essential difference with the results shown in 
Figs. 8(a) and 8(c) for the HgTe model. 
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Figure 8: (Color online) Localization length vs disorder 
strength U for the HgTe/CdTe model (a), (c) and for the Hal- 
dane model (b), (d). (a), (b) and (c), (d) correspond to ao 
and a x disorder respectively. In the calculations, the Fermi 
energy E is set to be E — 25meV, the topological mass is 
M — — lOmeV and the width of the ribbon is L y = 50a for 
the HgTe/CdTe model. The Fermi energy is E = 0.80i, the 
staggered term is A = — 1.60i and the width of the ribbon 
is L y = 40 x 3a for the Haldane model. The black line and 
green line are plotted for the ribbon and cylindrical geome- 
try respectively. Here, ao and a x correspond to the Anderson 
disorder and bond disorder respectively. 



IV. CONCLUSIONS 

In summary, we have studied the influence of bond dis- 
order (called a x disorder) on two models. This type of 



s 



disorder can originate in the deformation of the lattice, 
mismatch between two lattices or some chemical effects 
and thus exists extensively in many real materials. Un- 
like normal Anderson disorder, bond disorder does not 
induce a phase transition to a topological Anderson in- 
sulator (TA1). This is analytically verified by using Born 
approximation theory. The conductance and its fluctua- 
tions are then calculated, and the conductance plateau, 
which is the hallmark of TAI, is very good for Anderson 
disorder but becomes barely distinguishable for bond dis- 
order. In addition, phase diagrams are compared for two 
models with two types of disorder, Anderson disorder o~o 
and bond disorder a x . For Anderson disorder, the TAI 
phase can be seen clearly in the phase diagrams. How- 
ever, the TAI phase disappears for bond disorder. 

Because of a effective characterization in the metal- 
insulator phase transition, the localization length of the 
wave function is studied for the two models with Ander- 
son and bond disorders. For the two models with normal 
Anderson disorder, the localization length shows a peak 
at moderate disorder region, which exactly corresponds 
to TAI phase. However, the wave functions are localized 
quickly and no peak of localization length can be seen for 
the system with bond disorder. That means that Ander- 
son disorder and bond disorder can play different roles in 
the topological phase transition. 



To sum up, bond disorder can prohibit the system un- 
dergoing a phase transition to TAI, contrary to what 
may be seen in a Anderson disordered insulator.— These 
general conclusions in this paper are not restricted to 
two-dimensional systems and further work in a three- 
dimensional systems, e.g. B^Se-z model,— will be needed. 
If researchers intend to change the property of a topolog- 
ical insulator by addition of impurities, eg. shifting the 
energy level from the bulk states to the bulk gap, it will 
be necessary to guarantee that the impurities are mainly 
of the Anderson disorder type but keeps away from the 
bond disorder type. This is a key point for the doping 
process in preparation of topological insulators. 
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